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Small non-spherical perturbations of a spherically symmetric but time-dependent background 
spacetime can be used to model situations of astrophysical interest, for example the production of 
gravitational waves in a supernova explosion. We allow for perfect fluid matter with an arbitrary 
J?«.' I equation of state p = p{p, s), coupled to general relativity. Applying a general framework proposed 

by Gerlach and Sengupta, we obtain covariant field equations, in a 2-1-2 reduction of the spacetime, 
for the background and a complete set of gauge-invariant perturbations, and then scalarize them 
using the natural frame provided by the fluid. Building on previous work by Seidel, we identify 
a set of true perturbation degrees of freedom admitting free initial data for the axial and for the 
/ > 2 polar perturbations. The true degrees of freedom are evolved among themselves by a set of 
coupled wave and transport equations, while the remaining degrees of freedom can be obtained by 
quadratures. The polar i = 0, 1 perturbations are discussed in the same framework. They require 
gauge fixing and do not admit an unconstrained evolution scheme. 
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O ' I- INTRODUCTION 

a^■ 

^i^ \ In many situations of astrophysical interest, spherical symmetry is a good approximation for modeling a star in 

O I general relativity. One possible direction in which to go beyond that approximation is to allow for arbitrary linear 
perturbations, in order to add gravitational radiation to the picture. This introduces new physics, as the star can now 
5-H I lose energy through gravitational radiation, and a new window of observation, as this gravitational radiation can be 
P*J' detected. Gravitational wave detectors are expected to operate at the necessary sensitivity for the first time within 
^ , a few years, and a large effort is under way to model possible sources of gravitational radiation. If one allows the 
k>l spherical background solution to be time-dependent, as we shall do here, one can model for example the gravitational 

rS \ radiation emitted in a (slightly nonspherical) supernova explosion. 

jrt ' We assume here that the matter content is a perfect fluid described by an equation of state p = p{p, s) where 

■ ■ ' p is the pressure, p the total energy density, and s the entropy per particle. As a consequence of the perfect fluid 

approximation, s is assumed to be constant along particle trajectories, that is, we neglect the possible sources of 

entropy generation: heat fluxes, viscosity and chemical reactions. We also assume that there is only a single fluid 

present. 

Both the assumptions of approximate spherical symmetry and perfect fluid matter may be unrealistic for supernovae. 
Some supernovae are now conjectured to be quite nonspherical, and neutrino transport is believed to play an important 
role. Here we concentrate on giving a clean mathematical description of an almost spherical perfect fluid, in the belief 
that this approximation will be useful in some applications. 

There are many papers on the linear perturbations of a static spherical star, notably a series of papers by Thorne 
and coworkers ||l|-|6[|, another series by Cunningham, Price and Moncrief |7[], and a paper by Ipser and Price [^. 
The spherical symmetry allows one to decouple the perturbations into spherical harmonics. Because of the time- 
independence of the background, one can consider perturbation modes of the form exp(ia;t)/(r) and solve an ODE 
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eigenvalue problem for the mode functions f{r). The linear perturbations of a time- dependent spherical background 
were evolved in time by Seidel and coworkers ^. On a static background two of the perturbations obey trivial 
equations, and these have sometimes been overlooked in counting the degrees of freedom. The present paper sets up a 
mathematical and numerical framework for the study of arbitrary linear perturbations on a time- dependent spherically 
symmetric background spacetime. 

Even if one fixes the coordinate system (also called the gauge) in the background spacetime, the coordinate free- 
dom of general relativity poses a fresh problem when linear perturbations are added: one cannot easily distinguish 
infinitesimal physical perturbations of the background from infinitesimal coordinate transformations on the unchanged 
background. One can overcome this problem either by fixing the perturbation gauge, or by introducing linearly gauge- 
invariant perturbations. We choose the latter approach because it can easily be reduced to any particular gauge choice. 
The construction of the gauge-invariants, carried out in Sec. Ill, is a straightforward application of a general frame- 



work (for spherical backgrounds with arbitrary matter) due to Gerlach and Scngupta (from now on GS) |10| that is 
reviewed in Sec. II. 

The equations of GS are covariant in a natural 2+2 split of the spherically symmetric spacetime. These equations are 
only ready for numerical work after we have broken them up into evolution equations and constraints in Sec. IV. For 
the axial parity perturbations we can again use a general prescription due to GS. For the polar perturbations posing 
the initial value problem is difficult, and we use a key idea from Seidel ||9(] in order to find a subset of perturbations that 
can be given free initial data and that evolve among themselves, thus deserving the name true degrees of freedom. (The 
remaining perturbation components are obtained from these by solving the constraints.) We use the fact that the fluid 
matter provides us with a natural frame field (or set of observers) to decompose all tensors and tensor equations into 
scalars and scalar equations, and to distinguish evolution equations and constraints, without introducing coordinates. 

In many respects our framework is similar to that of Seidel and coworkers. However, our derivation of the gauge- 
invariants is more systematic, and we first make a covariant 2-1-2 split of the spacetime, then split the reduced 
2-dimensional tensors and tensor equations into frame components using the natural frame provided by the fluid. The 
final expressions are written so that they clearly display their causal structure. Being independent of background 
coordinates and linearly gauge-invariant, they can easily be specialized to a particular background coordinate system 
and perturbation gauge choice, while going in the other direction would be difficult. Seidel also restricts to ^ = 2 
angular dependence, and again it would be difficult to reconstruct the general from the special case. In the Appendix, 
we compare our notation and results to those of Thorne and Seidel. 

For completeness, the polar I = and I = 1 perturbations, for which the gauge-invariant treatment breaks down, 
are discussed here using equations that are as similar to the ones for Z > 2 as possible. In previous treatments, with 
the exception of Q , they were often neglected because they do not couple to gravitational radiation. 

Because the main application of this framework is to be stellar collapse, we must describe the matching of pertur- 
bations from the interior of the star (with both matter perturbations and gravitational waves) to the vacuum exterior 
(where the matter perturbations are meaningless). This will be done in a future paper. 

In the Appendix, we give the background field equations in polar-radial and comoving coordinates, discuss the 
special case of a static background, and compare our notation with that of Thorne and coworkers and that of Seidel. 
Finally, we describe a numerical algorithm for evolving the perturbations that we have successfully tested in an 



application to critical collapse |11 ij]. It is simple, second order, stable, and does not require special techniques at 
the center. As it explicitly uses the characteristic speeds, it allows us to impose purely outgoing boundary conditions 
at the outer boundary of the numerical domain. 

II. REVIEW OF THE GERLACH AND SENGUPTA FRAMEWORK 

A. The background spacetime 

In describing a spherically symmetric background spacetime and its linear perturbations we follow the route of 
reducing the system to 2-1-2 spacetime dimensions in a covariant manner, without introducing coordinates. In the 
following, we use abstract index notation, where Greek indices denote 4-dimensional spacetime, upper case Latin 
indices the 2-dimensional (H-l) reduced spacetime, and lower case Latin indices the orbits of the spherical symmetry 
(2-spheres). We write the general spherically symmetric spacetime as a manifold M = M^ x S'^ with metric 

5^^ = diag {gAB, r'^lab) , (1) 

where qab is an arbitrary Lorentzian metric on M^, r is a scalar on M^, with r ~ defining the boundary of IVP, and 
7at is the unit curvature metric on S*^ . Points in M^ are round spheres of area Anr'^ . (We can use r as a coordinate 



on M^, but do not have to.) Wc introduce covariant derivatives on spacetime, the subspace M^ and the unit sphere 
separately, with the notation 

9lii^;\ = 0, gAB\C = 0, Jab-.c = 0. (2) 

We shall also need the totally antisymmetric covariant unit tensors on M^ and S*^ respectively: 

— A BC _ B _ A be _ b /o\ 

^AB\C = 0, tAC^ = -gA , ^ab:c = 0, iac^ = 7a ■ (3) 

We parameterize the stress-energy tensor in spherical symmetry as 

tf_,^ = diag {tAB,Qr'^lab) ■ (4) 

With the shorthand 

VA = ^ (5) 

r 



the Einstein equations G^^ ~ 87rt^y in spherical symmetry, in the 2+2 split, are 

- 2{va\b + vavb) + (2wc''^ + ivcv^ - r-^)gAB - STrtAS, (6) 

wc''^ + wcw'^-7^ = 87^0, (7) 

where TZ = ^R a is the Gauss curvature of gAB- The equation of stress-energy conservation in spherical symmetry is 

tAB^^ + 2tABV^ ~2v'^Q = Q. (8) 

B. Nonspherical perturbations 

Any linear perturbation around spherical symmetry can be decomposed into scalar, vector or tensor fields on M^ 
times spherical harmonic scalar, vector or tensor fields on S*^. The spherical harmonic scalars on S"^ obey 

r'yr-.ab - -Ki + ^)yr- (9) 

Harmonic vector and tensor fields on S*^ can be constructed from the scalar harmonics. Wc shall need only vectors 
and symmetric tensors of rank two. A basis of harmonic vector fields on S"^ is formed for ? > 1 by 

A basis of harmonic symmetric rank-two tensors is formed by 

\^m r/ra \rra , \ ^ I \rra am , am /I 1 "\ 

^l lab, A afc=^i :ab^ ^ ^l lab, t>i a:b + ^l b:a^ (Hj 

where the last two expressions vanish identically for I = 0,1. Linear perturbations with different l,m decouple on 
a spherically symmetric background. In the following we consider one value of l,m at a, time, and no longer write 
these indices on Y, Sa and Zab- Furthermore, perturbations with different values of m for the same / have the same 
dynamics on a spherically symmetric background, so that m will never appear in the field equations. Perturbation 
fields containing an even power of eab, such as Y-a, are called polar, or even. They decouple from tensor perturbations 
containing an odd power of Sab, such as Sa, which are called axial, or odd. (Note that even and odd in this sense are 
not the same as even and odd parity in the standard sense.) 

The general axial metric and matter perturbations are parameterized as 



hT^'^Sa 

Symm h{Sa:b + St,:a) 



^3a"^— ( Q,,^^ h(Q . _i_ Q. ^ ) ' (-'^^) 



Symm At{Sa:b + Sb-.a) 
The general polar metric and matter perturbations are 



Ai^i,— c,^^^ A-ttQ , ^ Q, \ ] ■ ^^^> 



'^^''— [ Symm r^{KYjab + GY,,) ' ' ^'^> 



AtABY AtP°'"K„ 

Symm r'^At^Y'jab + AfZa 






We retain the notation of GS, except that we have added the superscript polar or axial where necessary to remove an 
ambiguity. Note also that r^ does not multiply AP in the last equation, and that At^ and AP are scalars. 

Let X be an arbitrary tensor field on the background spacetime, and AX its linear perturbation. Under an 
infinitesimal coordinate transformation x'^ -^ x^ +^'^, the perturbation is mixed with the background and transforms 
as 

AX ^ AX + C^X. (16) 

The perturbation AX is gauge- invariant to linear order if and only if £^X = 0. The existence of gauge-invariant 
perturbations is therefore linked to symmetries of the background solution. Because of the spherical symmetry of 
the background, the nonspherical perturbations with I > 2 can be made gauge- invariant, while I = and I ~ 1 
perturbations need to be considered separately. 

The general infinitesimal coordinate transformation can be parameterized as 

e^ EE (^UY, r^^Y.,a + rHiSa) . (17) 

GS work out how all the bare perturbations defined above transform. Then they form invariant linear combinations. 
Again we consider polar and axial perturbations separately. The split is as follows. There are 3 axial metric perturba- 
tions, and 1 axial infinitesimal coordinate transformation, leaving 2 axial gauge-invariant perturbations, in the form 
of a vector field on Af^. There are 7—3 = 4 polar gauge- invariant metric perturbations, in the form of a symmetric 
tensor and a scalar. There are 3 axial and 7 polar gauge-invariant matter perturbations. The axial gauge-invariant 
perturbations are 

l>l: kA^ hT""^ - h^A + '2hvA, (18) 

La = AtT""^ - QhT''\ (19) 

l>2: L = At-Qh, (20) 

and the polar gauge-invariant perturbations are 

? > : kAB = hAB - {PA\B + Pb\a), (21) 

k = K- 2v'^pA, (22) 

Tab = AtAB - tAB\cP^ - tAcP% ~ tBcP%, (23) 

t3 = At3 - (Qic + 2Qvc)p'' + ^-^^^QG, (24) 

1>1: Ta^ AtP°'" - tAcp"" - yQG\a, (25) 

l>2: T^ = At^ - r^QG, (26) 

where with GS we define the shorthand 

P^^<--VG|^. (27) 

These objects are not gauge-invariant for / — 0, 1, but it is useful to work with a single set of definitions and field 
equations for all I. We just have to take into account that for 1 — 0,1, G, h, T^ , L vanish and for / = 0, h^, h'^,TA, La 
also vanish, and we have to impose an additional gauge fixing in these cases. 

To obtain the bare perturbations in an arbitrary gauge, one fixes h, hF^ ^^ and G arbitrarily, and obtains all the 
other bare perturbations algebraically from these and the gauge-invariants. 

There is also a preferred gauge in which h = hF^ **' — G = Q, the Regge- Wheeler gauge |13 . The gauge- invariants 



have been defined to correspond one-to-one to the remaining bare perturbations in Regge- Wheeler gauge (from now on 
RW gauge). One can therefore describe the GS framework also as fixing the gauge to RW gauge. (The GS framework 



then still tells us how to transform to any other gauge.) The general axial metric and matter perturbations, in RW 
gauge but expressed through the gauge-invariants, are 

.RW_ / LASb 



^^^- - ( Symm L{Sa:b + S^-.a) ' ' ^^^^ 



The general polar metric and matter perturbations are, in RW gauge. 



A.gr- ( ^'f .J.. . ] , (30) 



kr^Y-ia 



b 

^V. - (^ Symm r^T^Yj^b + T^Z^b ) ' ^"^^^ 

These expressions are useful both in making calculations and in interpreting the results. 

III. THE PERFECT FLUID 

We now specialize the GS framework to a spherically symmetric perfect fluid coupled to gravity. In the perfect fluid 
approximation the entropy per particle s is constant along each particle trajectory (although it may vary between 
trajectories), and the pressure is isotropic. That is, we neglect heat fluxes, viscosity and chemical reactions, and 
assume that the fluid has only one component. The fluid is then completely characterized by a single two-parameter 
equation of state p = p{p,s). Within this class, we allow for arbitrary equations of state. Often it may be a good 
approximation to assume the entropy is constant in space and time, so that the equation of state is simply p ~ p{p) 
(barotropic fluid). This case can be recovered from our results by setting the specific entropy s to a constant and 
its perturbation a (introduced below) to zero throughout. A generalization of our framework to a fluid consisting 
of several (noninteracting) components would be straightforward: the relative abundances of particle types would be 
treated in the same way as s and a. 

Instead of giving p = p{p, s), many authors give p — p{n, s) and p — p{n, s) separately, where n is the conserved 
particle number density. The two specifications are completely equivalent. We have chosen the former as it is the 
simpler one from a spacetime point of view: n never appears in our equations. 

A. The perfect fluid background spacetime 

The perfect-fluid stress-energy tensor is, 

i^iv = ip + p)Uf^u^ + pgfj^^, (32) 

where p is the pressure, p is the density and u^ is the fluid 4- velocity. In spherical symmetry, u^ ~ {ua,0). This 
provides an orthonormal basis on M^, namely the timelikc unit vector ua and the spacelike unit vector 

UA = -eABU^- (33) 

We shall use this basis to transform all tensor fields and tensor equations on Af ^ into scalar fields and scalar equations. 
The metric and fundamental antisymmetric tensor on M^ can be written as 

gAB ^ -UAUB+nATLB, eAB^nAUB-UATlB- (34) 

For the spherically symmetric background the stress tensor becomes 

tAB ^ PUAUB +pnAnB, Q^P (35) 

We introduce the frame derivatives 

,/ - w^,/|A, ./' = n^f\A. (36) 



They will be useful later in discussing the initial value problem. In order to write our field equations using only scalar 
quantities, we introduce the following background scalars: 

n = \np, U = u^VA^r/r, W -n^VA^r' /r, h = u^\a, v = n^\A- (37) 

Note that fJ, ji and U are 0{r^) and even as functions of ?-, with U — ji = 0{r^), W is odd and 0{r^^), and v is odd 
and 0{r). 

It is useful to define the (spacetime-dependent) Hawking mass m in spherical symmetry as 

m ^ ^ (1 - r.Ar'') = '-[l + r\U^ - M^^)] (38) 

It is well known that the limit of m in a spacelike and future null direction gives the ADM and Bondi masses 
respectively, and that m = r /2 indicates an apparent horizon, m is odd and 0{r^) as one would expect. It can also 
be thought of as an integral over the density, in the sense that 

m' ^ATir'^{rW)p. (39) 

It is often useful to reparameterize U and W in terms of 

\v\^ = v\a = -U^ + W\ V=^- (40) 

Note that \V\ < 1, and in fact V is the velocity of the fluid with respect to constant r observers. |wp on the other 
hand is related to the Hawking mass via (|3^). 

The frame derivatives obey the following commutation relation: 

(/)'-(/')-=M/'-^/- (41) 

We also have 

UA\B ^nA{nBP~UBiy), nA\B=UA{nBi^-UBv). (42) 

The fluid equations of motion are given by the conservation of energy-momentum and by the assumption that the 
entropy per particle is constant along particle world lines. In our notation, the matter equations in spherical symmetry 
can be written as 

£7 + ( 1 + ^ j (2C/ + /i) = 0, (43) 

cln' + Cs +(l + -\v^ 0, (44) 

s = 0. (45) 

The first of these equations is the energy conservation equation (note that V^u'* = 2C/ + /i), while the second is the 
Euler, or force, equation. The third equation s = follows from conservation of the stress-energy tensor (p3) and the 
first law of thermodynamics. Here we have defined the shorthands 

2 _ (dp\ ^_ 1 (dp 



dp) ^ ' p \ds 



-i^ht > C^-[^ . (46) 



p 



We shall confirm later that c^ is the speed of sound by displaying a wave equation for sound waves. 
From the three Einstein equations (^) and the identity (^ we obtain the relations 

U'= Win - U) (47) 

W^ U{v - W) (48) 

W'^ -Airp -W^ + Up+— (49) 

[/= -Attp -U^ + Wiy~^ (50) 

From the identity g \ab — Oi using the fourth Einstein equation ([7|) to eliminate TZ, we obtain the useful relations 

p-iy' + fi^ -i'^ =n = ~4:Tr{p + p) + -^. (51) 

We now have a complete list of background identities that we can use later to simplify the perturbation equations. 
Expressions for the scalars in U, W, p,, v and m in specific coordinate systems are given in the Appendix. Note that 
[/ and p vanish on a static background. 



B. Nonspherical perfect fluid perturbations 

There should be 5 independent fluid perturbations, namely a density perturbation, entropy perturbation, and a 
3-velocity perturbation. It turns out that 4 of these are polar and 1 is axial. The polar perturbation of the fluid 
4-velocity is 



while the axial perturbation is just 



iriA + T:hABu" ) Y, aYa 



AUf, = {OjSa). 



(52) 



(53) 



Note that the ansatz for Au^ ensures that A{g^^u^^u'^) = to linear order, a and (3 parameterize polar and axial 
tangential fluid motion, while 7 parameterizes perturbations of the radial fluid motion. The density, pressure, and 
entropy perturbations are 



Ap = LopY, As = o-y, 

A gauge-invariant set of fluid perturbations is 



Ap = ^Ap + ^As = (c^c^ + Ca)pY 
op OS 



a = a — p ub, f3 = 13, 7 = 7 — 71 



P^ua\b + 7;'^^{PB\A - Va\b) 



P^^\A, 



(7 — a — p siA 



(54) 

(55) 
(56) 



with PA as defined above. 

For the gauge-fixed viewpoint, we give the fluid perturbations in RW gauge, expressed through the gauge-invariant 
perturbations. The polar perturbations are 



AwJ^w 



illT-A + -kABU^)Y, aYa 



Ap^^ = wpF, 



A, 



,RW 



== aY, Ap«w = {ciu; + Ca)pY, 



and the axial perturbation is 



Au 



f^ = iO,PSa). 



(57) 



(58) 



The general gauge-invariant stress-energy perturbations of GS can be expressed in terms of the gauge-invariant fluid 
perturbations we have just defined. In the axial sector, we have 



La = f3{p+p)uA, 
and in the polar sector, 

Ta = ct{p + p)uA, T^ ~ pk + clpuj + Cpa 

Tab = (p + p) 



L = 0, 



(59) 



1 



jiuAriB + TIaUb) + ^{kACUB + UAkBc)u^ 



T^ = 0, (60) 

+ Ujp{uAUB + cIuaUb) + aCpUAUB +pkAB- (61) 



The only axial fluid perturbation, /3, describes equatorial fluid rotation, uj describes total density perturbations, a 
entropy perturbations at constant total density, 7 describes radial fluid velocity perturbations, and a describes the 
velocity of tangential fluid motion between the poles and the equator (azimuthal displacement). 



IV. THE INITIAL VALUE PROBLEM FOR THE FLUID PERTURBATIONS 



The previous section has been a straightforward application of the GS formalism. The perturbed Einstein equations 
and stress-energy conservation equations are given in terms of the gauge-invariants and in 2-1-2 covariant notation by 
GS, and we only need to substitute our parameterization of the gauge-invariant stress-energy perturbations into their 
equations. Note that for the perfect fluid the conservation of stress- energy, plus the conservation of entropy along 
particle world lines, is equivalent to the equations of motion. (For a barotropic fluid with p = p{p), stress-energy 
conservation alone is sufficient.) We now turn to the less trivial aspect of breaking these equations up into evolution 
equations and constraints, that is, of posing and physically interpreting an initial value problem. For the axial sector 
this is an application of a general method described by GS. For the polar sector we use an idea due to Seidel H to 
isolate the true degrees of freedom. 



A. Z > 2 axial perturbations 



Extracting a well-posed initial value problem from the gauge-invariant perturbation equations for the axial sector 
is relatively straightforward. The one nontrivial matter conservation equation is (GS4) 



{r'LAy^^il-m + 2)L, 



(62) 



Note that for perfect fluid matter, we have L = 0. In our notation, and using (|43 ) to eliminate ft, the (3 equation of 
motion is 

P ~ cl{2U + fi)(3 ^ 0. (63) 

If the perturbed fluid velocity is to be a regular vector fleld at the origin, [3 must scale as /3 ~ r'+^. We define 

(3 = r'+i/3. (64) 

The new variable P can be expanded in even positive powers of r near r = 0. Its equation of motion is 

P + [-cl{2U + fi) + {I + 1)U] P = 0. (65) 

As we have L = 0, Eq. (|2|) can be integrated to 

r^A = eABT^"", (66) 

where T is a scalar which contains all the information about the axial matter perturbation. The two components of 
this equation are 



T' = -r2(p+p)/3, r = 0. 



(67) 



[The integrability condition for these equations is just (|6|), using (pl|).] That is, we set /3 on a given time slice and 
then we integrate the T' equation for T on that slice. The T = equation states that T is constant along the integral 
curves of u"^, that is, each particle of the fluid sees a constant value of T around it. In this sense, the axial matter 
perturbation is time-independent, even on a time-dependent background. 
GS have shown that by defining the scalar 



n^e^^{r-'kA) 



(68) 



the Einstein equations for kA can be reduced, for any matter content, to the scalar wave equation valid for I > 2 
(GS17, corrected GSII6.5b') 



-{r'n)\A 



{I - l){l + 2)n = -16TTe^^LA\ 



(69) 



which they call the odd-parity master equation. In the case of perfect fiuid matter, the equation of motion of La does 
not contain 11, and so can be solved on its own. 

The behavior of fc^ at r = can be calculated by going to Cartesian coordinates. In order to enforce the correct 
scaling of kA, we define 



n EE r'-^n. 



(70) 



where the rescaled field 11 can be expanded in positive even powers of r at r = 0. For 11 we obtain the field equation 



n 



2(1 + l)v^niA + (l + 2) 



r\A 



\A 



2m 



a-1)- 



n 



-16nr{p + p){ /3' + 



(/ + 1)W - 



n' 



i + ? 



Near the origin, this equation has the approximate form 

- fl + U" + 2{l + l)r-^U' + exphcitly regular = 0, 



/3 } ■ (71) 



(72) 



which is suitable for numerical solution. H' is an odd function of r, so that H'/r is regular at r = 0. 

The scalar 11 codifies all the information about the evolution of the axial parity gravitational wave. In fact, once 
we know 11 we can reconstruct fc^ using Eq. (GS16c): 



{l--l){l + 2)kA--16TTr^LA-eAB{r^ny^ ^eAB{l67rT-r^ny , (73) 

where the last equation holds because of L = 0. This equation is valid for I > 1. (For I = there are no axial 
perturbations.) For I > 2, we obtain kA from 11. For / = 1, however, we obtain a constraint on 11 instead. The case 
I — 1 has to be treated separately. 

The free axial perturbations (for / > 2) are the fluid velocity perturbation /3, which obeys a transport equation 
independent of 11, and the metric perturbation H, which obeys a wave equation with /3 as a source. On an initial 
Cauchy surface, one can freely specify /3, 11 and 11, and in this sense there are three (first-order) degrees of freedom. 
Note that even though /3 itself cannot oscillate, a non- vanishing /3 can couple to nonlinear oscillations of the spherical 
background to generate waves in 11. 

Ji. I — 1 axial perturbations 

Physically, we cannot have any degree of freedom in the metric for I — 1 and therefore 11 cannot obey a wave 
equation. It has to be constrained. This is clear in Eq. d73), which for / = 1 implies: 



r'^U = 16ttT + const. (74) 

The integration constant we have introduced here parameterizes a Kerr-like angular-momentum perturbation that is 
singular at the origin. It is zero if we demand a regular center. Note that even for Z = 1, 11 is still gauge-invariant, 
but kA cannot be reconstructed uniquely form it. Now we have to invert Eq. (pq), obtaining kA with a gauge freedom 
that is the gradient of an arbitrary scalar. 

C. I > 2 polar perturbations 

The polar perturbations are more entangled than the axial ones. Here we follow Seidel B in first focusing attention 
on those components of the linear Einstein equations that do not contain matter variables (except for the entropy 
perturbation). We shall see that these Einstein equations (plus the evolution equation for the entropy perturbation) 
can be solved autonomously: derivatives of the metric perturbations contain all the information about the matter 
perturbations (except for the entropy perturbation). 

We use the fact that the fluid provides a natural frame in Af^ to decompose the symmetric tensor kAB m a 
coordinate-independent way into three scalars, as 

kAB = Vi-UAUB + nAUB) + (J){UAUB + UaUb) + IpiuATlB + TIaUb)- (75) 

The seven linearized Einstein equations of polar parity can be taken as (GSlOb), (GSlOc), (GSlOd) and (GS12, 
corrected GSII3.13a). We go over to a system of scalar equations by substituting the definitions (JT^), (|3|) and ( |37| ) 
into terms like kABiCi ^'^d projecting the result on the basis (u, n). The left-hand sides of the perturbed Einstein 
equations decomposed in this way are given in p2[ (with respect to an arbitrary frame). Here we give only the 
components of the right-hand sides for the perfect fluid (with respect to the fluid frame): 

l>0: u^n^TAB = -{p + ph+^{p-pH, (76) 

M^M^TAi3=pc^ + p(r;-0), (77) 

n^'n^'TAB = p(c?w + Ca) + p(r; + </.), (78) 

T^ :^ pic^uj + Ca) + pk, (79) 

1>1: u^TA = -{p + p)a, (80) 

u^Ta = 0, (81) 

l>2: T^ = 0. (82) 

We have indicated for what values of I the corresponding Einstein equations are valid. Because u is always timelike 
and n"^ always spacelike the scalar decomposition allows us to disentangle evolution equations and constraints. 



In the remainder of this subsection we discuss the case I > 2, but we always indicate which equations are also valid 
for / = 1 and / = 0. The dynamics oi I = 1 and I = are discussed in the following two subsections. 

For I > 2 the vanishing of the source T^ implies r] = 0. Among the other six source terms we have just given, there 
are three linear combinations that do not contain a, to, or 7, and these give rise to three perturbed Einstein equations 
not coupled to those three matter perturbations (they do contain the entropy perturbation a) . If we choose the three 
combinations to be 

T^ + 2{n^TAy - n'^n^TAB + 2(2i/ - W)n'^TA, (83) 

{-clu^u^ + n'^n^)TAB + 4VFn^rA, (84) 

n^TA, (85) 

and introduce the new variable 

X = <f) — k-\-rj, (86) 

which from now on replaces (j) as an independent variable, then the corresponding Einstein equations take the form 

1>1: -x + x" + 2(M-t/)^'-^x. (87) 

~k + clk" -2clUij'^Sk, (88) 

-i^^S^. (89) 

These are the core equations of our paper. Only the highest derivatives of all variables, which determine the causal 
structure of the equations, have been written out explicitly. The source terms S^ , Sk and 5'^ are homogeneously linear 
in X and k, and their first derivatives x', fc', %, and k, as well as ip and a (undifferentiated), with background-dependent 
coefficients, but do not contain any higher derivatives. They are given in Appendix H. 
We choose the remaining three components of the Einstein equations as 

(90) 
(91) 
(92) 

which give us 

;>0: 8TT{p + p)j^{ky + C^, (93) 

a, (94) 

1>1: 16Tr{p + p)a ^ ij' + Ca, (95) 

where again all the highest derivatives have been written out explicitly, and the right-hand sides do not contain a, 
7 or w, (or indeed a). These equations give a, 7 and uj as spatial derivatives of the metric perturbations. Seen the 
other way around, they form an ODE system for i/j, k and k in the radial coordinate which can be solved from the 
center outwards, with boundary conditions at the center given by local flatness. The correspondence between a, 7 
and (jj on the one hand, and ip, k and k on the other (for given data x and x, and a) is therefore one-to-one. (With 
one proviso discussed at the end of this section.) 

The linearized stress-energy conservation equation (GS15a,b, corrected GSII3.14a,b) are equivalent to the following 
equations of motion for the matter perturbations. 

l>0: -u;-(l + ^) G+ 2) "^"' ^^^^ 

^)(7-|)-+cV = 5„ (97) 

1>1: -d^S'a. (98) 

Again we have written out only the highest derivatives. The three source terms are linear in the (undifferentiated) 
matter perturbations a, 7 and to, and a, as well as %, k, V 1 k' , x, k, and tjj. The full equations are given in Appendix 
^. It is possible to remove tp and tp' from (96) and ( |97| ) using the Einstein equations ( |89| ) and (|95|). The matter 



(u^n^ 


''+n\ 


'')Tab, 




u^ 


u^Tab, 
u^Ta, 


Snip 


+ P)l ^ 


--{ky + c„ 




Snpuj - 


- -k" + 2Uij' 


L67r(p 


+ p)a = 


--^' + c^, 



perturbation equations then do not contain highest derivatives of the metric perturbations. In this form they are 
most amenable to numerical solution. We obtain 
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/>1: -d,-(l + ^U' = S^, (99) 

l + -J7 + cy = 5^. (100) 
To these matter equations we have to add the perturbation of the trivial entropy equation s = 0, which is 

l>0: (7+ (7+1 j .s' = 0. (101) 

We now introduce auxihary variables 'x' ,k' ,x^k in order to consider the perturbation equations as a system of 
first- order system of evolution equations and constraints for the variables 

u = {x, ip, k, a, 7, uj, a, x' , k' , x, k}. (102) 

We consider an equation that contains ii an evolution equation, and an equation that contains only u and u' a 
constraint. Constraints can be solved within a single timelike hypersurface. (If u'^ is not normal to the surface of our 
choice, we need to use the evolution equations in order to write the constraints as ODEs in the radial coordinate.) 
Apart from the equations we have already given explicitly, the first-order form introduces trivial evolution equations 
for X s-nd fc, and trivial constraints for x' ^-nd k' . Altogether we have 11 evolution equations and 5 constraints for 
11 first-order variables u. This means that there are, in a first-order in time sense, 6 true degrees of freedom, or 6 
functions of the radial coordinate that can be specified freely on a Cauchy surface. A natural choice of these would 
be the matter perturbations a, 7, w and a, plus x and x which describe polar gravitational waves. [There are two 
"polarizations" of gravitational waves, namely axial (II) and polar (x).] The metric perturbations tp, k and k would 
then be obtained from the constraints. While this scheme is the most natural in terms of a split into matter motion 
and gravitational waves, it maximizes the number of constraints that must be solved. This is an inconvenience both 
in terms of numerical work and numerical stability. 

Seidel H has noticed that one can identify the 6 true degrees of freedom with the initial data for the metric 
perturbations alone, plus the entropy perturbation a. This means that 

Ubcc= {x,i^,k,(T,x,k} (103) 

can be set freely and evolve among themselves. In a first-order formulation, there remain the trivial two constraints 
and two evolution equations for the auxiliary variables x'j k' , x and fc, but the three nontrivial constraints are now 
only used if one wants to reconstruct a, 7 and u as spatial derivatives of the free Cauchy data Wfioc- 

Let us now consider the causal structure of the equations from a spacetime point of view. The highest derivatives of 
X form a wave equation with characteristics given by the metric g^^, and as we have seen, initial data x and x can be 
set independently from the matter perturbations. Therefore, x can reasonably be said to parameterize gravitational 
waves inside the matter. The two first-order evolution equations for to and 7 are equivalent to either of the two 
second-order equations 

-7 + ch" + --0, (104) 

-w + c^cj" + ...= 0, (105) 

where we have only written out the highest derivatives. These are wave equations with characteristics given by the 
"fluid metric" 

9^7'^ = diag {-cluAUB + riAriB^r'^jab) (106) 

The characteristics have speed Cg relative to the fluid background. The evolution equation for fc is a wave equation 
with the same characteristics. It "knows about" the speed of sound because c^ appears in the linear combination, Eq. 
(p4), of Einstein equations that gives rise to this equation. It is surprising that the metric perturbation k describes 
sound waves. The analogy is correct, nevertheless, as the data w and 7 are essentially the spatial derivatives of the 
data fc and fc. It may be helpful to think of fc as describing "longitudinal gravitational waves" , made physical by the 
presence of matter. Finally, "0 is advected with the fluid, as its evolution equation contains ip but not "0'. The same 
is true for the matter perturbation a, and not surprisingly, a is essentially the spatial derivative of ip. a describes 
tangential fluid motion between the poles and the equator, and is therefore coupled to the density perturbation w 
that acts as a restoring force. 
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Finally, we consider the behavior of the polar perturbations at the origin. We require that the metric perturbation 
( pOl ) be a regular tensor in four dimensions at r = 0, as defined by going to Cartesian, Minkowski-like coordinates 
near r — 0. We find that any regular solution must scale at r = as 



I > 



X- 



r'+^X, 



tP 



.'+1 



^, 



r'/c, 



V 



r'?7. 



(107) 



where the barred variables can all be expanded in even positive powers of r. This behavior at the origin is consistent 
with the equations of motion, in the sense that Eqs. (p7|), ( |8q ) and ( p9| ) can be solved for x, k, and -ip order by order 

in r, for arbitrary values of (%, x, k, k, ip, a). 

To investigate the regularity of the velocity perturbations at the origin, we work in RW gauge, and introduce 
Minkowski-like coordinates t,x,y,z near the center. For simplicity of presentation we only consider the case m = 0. 
The four-velocity perturbation becomes 



Au^'^dx 



^-2 



Y, 
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dt 



1 dYi^ 



7- 



-'Y 



m 



ar 



sin 61 de 
_2Cob6 dYio 



sin 6* d0 



dz 



[xdx + ydy + zdz) 



(108) 



\l-2^ 



This is regular if and only if the terms in square brackets are regular scalars. Now cot 6 dYio/dO = — ^y;o + 0[(cos ( 
and so the leading orders of a and 7 can both be a factor of r^ lower than naively expected if they are correlated. 
Taking this possibility into account, the matter perturbations near the origin behave as 



l> 1 



7 



„i-l; 



7, 



with 7 = ;« + O(r^) at r = 0, 



LU = r u), 



(109) 



To avoid the constraint between the leading orders of 7 and a, one could replace 7 by 7— a' as a dependent variable that 
would be unconstrained, similar to the variable x- The present notation, however, gives rise to the more transparent 
evolution equations. Again, the scaling at the origin is consistent with the equations of motion (|9^), (|9^) and (p^). 
In particular, what looks like a leading term that is too large by a factor of r~^ cancels in each case. Furthermore, 
the constraint between the leading terms of a and 7 is conserved by the evolution. 

The constraint equations (B3), (P4), (pa) can be solved consistently for 7, u) and a, order by order in r. The converse 

is almost true. For given x, x, a, a, tu and 7, these equations can be solved order by order for ip, k and k, with 
the exception of the leading order ko of k and the leading order i/Sq of i/j. The gravitational wave degrees of freedom 
are therefore the functions Xj X, plus the two numbers ko and ipo, while the matter degrees of freedom are either the 
functions to, 7, a, or alternatively the functions k without fco, k, and ip without tpQ. In practice there is no need to 
distinguish between matter and metric perturbations. 



D. Z = 1 polar perturbations 



The case of / = 1 for the polar perturbations differs from / > 2 in two important aspects. First, the trivial field 
equation that stated that 77 = no longer holds, so we have an additional variable. Second, the variables we have 
defined for I > 2 exist for / = 1, but are only partially gauge- invariant. Under a gauge transformation generated by 
the vector £,^dx^ — ^AYdx"^ + r^^Y-adx"^, we have the following change in the metric perturbations: 



kAB > kAB + (r^^\A)\B + {r^£,\B)\A, 

k^k + 2^+{r^)^%A- 



Decomposed in the fluid basis, this is: 



7] 

X 
k 



1] + 2r^{-Ui + WO + r2(-e + e" + v£.' - d), 



X - 2e + 2r' 



WS,'+£," + {U-lj)i 



^ ~ 2r'{Ue + Wi) - r2 [(0' + (C'J- M?' " ^C 
k + 2^ + 2r^-Ui + WC). 



(110) 
(111) 

(112) 
(113) 

(114) 
(115) 



We need to impose one gauge condition in order to fix the gauge freedom parameterized by ^. It is tempting to e xtend 
Seidel's free evolution scheme to I — 1. Clearly, one can use the gauge freedom to set 77 = by solving (112) as a 
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wave equation for ^. This still leaves us with two free functions to specify the gauge completely, such as ^ and ^ on 
one spacelike slice. As there are no dipole (/ = 1) gravitational waves, one would also like to make the gauge choice 
X = 0, as we have seen that x obeys a wave equation at the speed of light. While the two free functions could be 
used to set x and x to zero on one slice, x could not be made to disappear on every slice. Conversely, one might be 



able to solve (113) as a second-order constraint equation for ^ on each slice in order set x = 0; which would leave us 
with one free function of one variable, such as t] at the center as a function of time, {rj' at the center would be fixed 
by imposing regularity). That means that we would be stuck with at least this constraint to solve. Furthermore, the 



presence of the term ([/ — /x)^ in (113) complicates the interpretation as a constraint, as the hypersurface on which 
this equation becomes an ODE in the radial coordinate is not a natural one to choose for any other purpose. 

All this is not promising. Instead we follow Thorne Q in choosing the gauge fc = 0. We might call this gauge the 
radial perturbation gauge, as th e ba ckground scalar r remains the area radius in the perturbed spacctime. We now 
use the matter equations (|9q)- (|l0l| ) to evolve the matter perturbations w, 7, a and a. These evolution equations 
contain 77, ip and x (but not their derivatives), while k vanishes by our gauge choice. We can solve the five Einstein 



equations (88), (|89|), (|93|), (M), (|95[) for x and x', "0 and ■0', but only for one linear combination of f] and 77', namely 



Drj, where the differential operator D is defined as 



r^^fl 



Df = r-^\v\-^{Wr-Uf) = --^ (116) 

which is just d/dr on a spacelike surface that is everywhere orthogonal to the surfaces r = const, (a polar slice). 
The expression for Drj is given in Appendix |A|. The remaining sixth Einstein equation ( p7|) does not tell us anything 
new: as it contains x, it follows as a consequence of the other equations. The fact that we do not have r] and 77' 
independently means that we can integrate the constraint for 77 only on a surface that is normal everywhere to the 
surfaces r = const. As we are restricted to these surfaces anyway, we also give only Dx and Dtp in Appendix |A[ 

The three constraint equations form three coupled first-order linear ODEs. To analyze them we define again 
regularized variables: 

X = r^X, tp = r'^-ip, T] = rfj, (117) 

7 = 7, a = ra, with 7 = a -|- 0{r ) at r = 0, w = rui, a = ra. (118) 



The constraints given in Appendix |A| give two relations among the leading terms of these new variables: 

-2U'^fi + 0{r'^) atr = 0, (119) 



Wa\l + - 



tp = 8Tra{p + p) + 2Uf] + 0{r^) at r = 0. (120) 



OJ 

V P 



Note that imposing fc = we still have the residual freedom of functions ^ obeying ^ -t- r^\v\'^D£^ = 0. We have to 
give a boundary condition for this ODE at r = at each instant of time, and therefore the residual freedom is a single 
function of time. We can use this function to set fj = 0{r^) at the center. 

E. Z = polar (spherical) perturbations 

The case of spherical {I ~ polar) perturbations has the same problem with 77 as the I — 1 case. It differs in that 
the polar matter perturbation parameterized by a does not exist, and that there are now two (polar) gauge degrees 
of freedom, with no remnants of gauge-invariance. We impose the gauge 



2UW 



0' ^ ^^- 77^-1772 (^-x)-0. (121) 



Here ip \s ip in the radial (instead of the fluid) frame, see Appendix ^. In this gauge we have two metric perturbations, 
77 and Xi which obey constraints given in Appendix H. As in the I = 1 case, only the combination Dij of 77 and rj' 
is known. Using these constraints, the metric perturbations can be calculated from the matter perturbations 7,a;,o'. 
This gauge is the only one in which we have found simple evolution equations. Note that a natural coordinate choice 
for the background spacetime are polar-radial coordinates which are defined by ggg — r^, gtr = 0. I? is then just d/dr 
in these coordinates, and the gauge choice k — tp = just means that the perturbed metric (which is still spherical) 
has the same form as the background metric. Here we give the spherical perturbation equations in our notation, 
rather than in polar-radial coordinates, for the purpose of a unified presentation of all perturbations. Given that the 
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spherical perturbations are so much messier than the / > 2 perturbations, however, one might as well calculate them 
in the same gauge choice one has already adopted for the background spac etime. 



The matter perturbations obey the evolution equations (puj), (|97|), (101), which contain derivatives of the metric 
perturbations. The latter can be eliminated, using the perturbed Einstein equations. The resulting equations are 
given in Appendix ^. 

Again, we introduce regularized variables: 

X = r^X, 1 = ^, (122) 

7 = r7, to = uj, a = a. (123) 

Note that now we do not have any cancellation of leading orders in r, so we obtain the powers of r naively expected. 
The constraints give only one condition at the center: 

X = ^^ui - 2U^fi + 0{r^) atr^O. (124) 

Again, we have a residual gauge freedom that wc can use to set fj = 0{r^) at the center. 

V. CONCLUSIONS 

We have given the field equations for all spherical and non-spherical linear perturbations of a spherically symmetric 
but time-dependent self-gravitating perfect fluid. In this task we have applied a formalism created by Gerlach and 
Sengupta [fo|. We have generalized and clarified previous results by Thorne and coworkers ||l|-§[, Ipser and Price ||] 
and Seidel iP. Our formulation is distinguished by the following nice features. 

1) Our perturbation variables are linearly gauge-independent (with the exception of the polar I — Q,! perturbations, 
see below) . This means that we can be sure to only count physical degrees of freedom. Our results can be translated 
into any particular gauge (for example Regge- Wheeler gauge) by explicit algebraic formulas. 

2) Using the fact that the fluid provides a preferred frame, we have decomposed both tensor variables and tensor 
equations into frame components. This procedure clarifies how many degrees of freedom are present, which perturbed 
Einstein equations are evolution equations and which are constraints, and what are the characteristic speeds. We 
find that in the sense of free functions on a Cauchy surface, the axial I = 1 perturbations have 1 gauge-invariant 
degree of freedom, the axial I > 2 perturbations have 3 (including one gravitational wave polarization), the polar 
I = perturbations 3, the polar I = 1 perturbations 4, and the polar I > 2 perturbations 6 (including the other 
gravitational wave polarization and sound waves). 

3) All the final background and perturbation equations are given in terms of frame derivatives of spacetime scalars 
or frame components of tensors. This means that our equations can be easily translated into any coordinate system 
for the background (for example comoving or polar-radial coordinates). 

4) With the exception of the polar / = 0, 1 perturbations (see below), we have split the gauge-invariant perturbations 
into a set of variables that can be specified freely on a Cauchy surface, and that evolve freely among themselves, and 
another set of variables that are determined by the first set through algebraic constraints. Writing out only the highest 
derivatives of each variable, the axial perturbation free evolution equations for I > 2 are 

- ri + n"= . . . , 

/3= . . . , 

and the polar free evolution equations for I > 2 are 

-X + x" + 2i^l-U)^P'=..., 
-k + clk" -2clUil/:= ..., 

v>=..., 

where c^ is the local sound speed, and the coefficients U and /i vanish in a static background. Having a set of free 
evolution equations that are manifestly hyperbolic is helpful for numerical work, in particular if one wants to use 
characteristic methods. (An example is given in the Appendix.) 

5) No gauge-invariant variables exist for the polar / = 0, 1 perturbations. Nevertheless we define the same variables 
for Z = 0, 1 as for I > 2. As these variables are not fully gauge-invariant for / = 0, 1, additional constraints must be 
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imposed to fix the gauge, but we can make use of the fact that our variables are still partially gauge-invariant for 
1 = 1. This clarifies the role of the polar I = 1 perturbations, and in particular the number of degrees of freedom. 

6) We allow for a two-parameter equation of state p = p{p, s) and take into account perturbations of the entropy s. 
Our framework can easily be generalized to multi-component fluids. 

We believe that our formulation characterizes the perturbation initial value problem more clearly and concisely 
than previous work, and is also a good starting point for its numerical solution. 
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APPENDIX A: FULL POLAR PERTURBATION EQUATIONS 



Here are the missing source terms in Eqs. (|87|)-(p9|) and (|93|)-(100). Recall that some of the equations to which 
these sources belong are not valid for / = or Z = 1. On the other hand, note that rj = for / > 2. Note also that we 
have used the background equations to eliminate all occurrences of tl, Jl', [/, [/', W, W, and /i — i^' . 



2z/2 + Sttp - ^ - 2[/(/i - U) 



(, + ^ + (i^i)([±S, 



+3nx + 4(/i - U)k - {5i^ - 2W)x' - i^Ifiv - 2(/i - U)W + ^i' - i>]^p 

2 ^(^ + l) + 8 



+2?7" - 2{fi - U)fi + {8iy - 6W)7^' - 



-Au' 



SvW + 4:{2fiU + U^ - AW^ 



Sk = (1 + ci)Ux + m + clifi + 2U)]k - 1^(1 - cDx' -{i^ + 2Wcl)k' 



2 ^-W 



- Snp — c^ 



l{l + l) 



+ 2U{2^ + U)- Slip 



(x + fc) 



(/-!)(/ + 2) 
2r2 

^2C/r) + 2Wr]' - 



(1 + cl)x + 2[-pW{l - cl) + {v + W)U{1 + cl)]^ + SnCpa 



l{l + l) 



6W^ + Wnp - 2C/(2/i + U)cl 
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-X + 2[i^U+ip + U)W]ilj 



+Ux + (m + 2C/)A: 4- Wx' - 2Wk' - 2f]U{2fi + U), 
Ca = 2p{x + k) + 2i^^P + x + 2k- 2r]{p + U), 
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(Al) 



(A2) 
(A3) 
(A4) 



(A5) 
(A6) 



(A7) 

(A8) 
(A9) 
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5.. = ( 1 + - 

p 



1(1 + 1) , ^^ k , , , ,^ 

- 2 + 8^(/9 + p) a + ;r + (a* + C/)f? - m(x + fc) 



+ (a. + 2C/) ( c2 - ^ ) u; 
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(All) 



The matter perturbation equations for the polar I == 1 perturbations are the same as for ^ > 2, namely (|99[), ( 101 ) 
The three constraint equations for the polar I = 1 perturbations that we require are 



r\v\ Drj = Anpil ^ Cg)u} — 167r(p + p)Ua — AupCa 

- i^-iW'^ + ^TTp + uAr]- [W^ + U^ - 4:7r{p + p)] X - 2UWip, 

2 



r\v\ Z^x = STrpw — 327r(p + p)C/a - 



+ 2[/^ - Svrp X + 2 [t/i/ -{fi + U)W] iJj - 2U^r], 



r\v\^Dij = 87r(p + p)(7 + 2Wa) + 2{pW - vU){r] - x) + -WW-q - 



C/^ - W^2 + ^ + 47r(p - p) 



v-. 



(A12) 
(A13) 
(A14) 



The matter evolution equations for the polar Z = 0, or spherical perturbations are 



w- ( 1 + - )7' = (/i + 2[/)C(T - w 
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(A15) 
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The remaining two metric perturbations rj and x can be obtained from the constraints 



, ,2^ . / M 2[/2 \ ^ ^ ^2t/VK ^ C/2 + 1^2 

r\v\ Dtj = 4tt(p + p) [ x+ T-p'^ ) + 87r(p + p)-p-|2-7 + 47rp- 



2 (Ca+(l + c2H, 



(P + P) 
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APPENDIX B: BACKGROUND EXPRESSIONS IN POLAR-RADIAL COORDINATES 



Excepting null coordinates, there are two natural coordinate systems for spherically symmetric fluid spacetimes. 
Polar radial coordinates r, t use the scalar r as one coordinate, and make t orthogonal to it, so that the metric is 



ds' 



-a^df + a^dr"^ + r^dfl^ 



The gauge is fixed completely when a is fixed to be 1 at r = (or at r = cx)). In these coordinates we have 

, ,o 1 2m 



r 



r 



,t 



1 



V 



u'' = 



V 
1 



a 



VI -1^2' 



fl = U 



at . VVt 



V = u 



a 1 - y2 



fVa^r Vr 



\ a 1-v^ r 



a 



VVr 



The first derivatives of the metric that appear here are given by the Einstein equations 

a^ - 1 a, r ^ p + V'^p 

= Aira r- 



2r a 

1 — a^ a,r 



2r 



l-y2 ' 

— 47ra r —, 

a 1 - y2 ' 

a,t .2 V{p + p) 
— 4:TTa r T^. 

a i-y2 



(Bl) 

(B2) 
(B3) 
(B4) 
(B5) 
(B6) 

(B7) 
(B8) 
(B9) 



The Einstein equation ([7|) is a combination of derivatives of these three. The matter evolution equations are, elimi- 
nating all derivatives of the metric variables: 



a a a 



2V 
Aira^rip + pfV - {p + p){ — + Vr) + VCps ^(1 - V"^) 

r 



a a a 



^^na\{pclV^ +p)-{l~ l^2)P^^V + c|p^ ^ l_^ ^ ^-±^clV^ 



p + p 



2r 



2r 






0. 



(BIO) 

, (BII) 
(BI2) 



APPENDIX C: BACKGROUND EXPRESSIONS IN COMOVING COORDINATES, AND COMPARISON 

WITH SEIDEL 



The metric in comoving coordinates, using the notation of Seidel H, is 

ds^ = -N^dT^ + A^dp^ + R^dVi^. 



(CI) 



Note there are three independent metric functions, and that Seidel's R is our r. The condition that lines of constant 
/i are particle worldlines leads to 
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u^^AT-i, u'^^O, n^ = 0, nf'^A-^. (C2) 

This, and fixing the value of N at the center fixes the coordinate freedom, up to a relabeling of the radial coordinate 
^ — > fJ-'{fJ-). In these coordinates we have (using our notation on the left-hand side, and Seidel's on the right-hand 
side) 

^=A^' ^=Ai?' ^^AAf' ^=AiV' - = ^-^ + ^- (C3) 

Seidel fixes the gauge freedom in /i by imposing A = {AiTpR^)^^. This means that his radial coordinate is mass-like, 
with /i ~ Anr^ p/3 to leading order. 

The gauge-invariant metric perturbations of Seidel [p[ are related to our variables in the following way. (Seidel's 
notation is used on the left-hand sides, and our notation on the right-hand sides.) 

qi = + 7?, (C4) 

92 = fc, (C5) 

N-^TTi = (0 + f])-+2ilj' + 2v%l^ + ^(0 - r/), (C6) 

N-'^'K2 = k + 2Wij + U{(j)-T]), (C7) 

iV-M-ii?2^4 = tj}. (C8) 

The r; terms are contained in Seidel's definition, but 77 = for / > 2 by one of the Einstein equations. Seidel's 
conditions for re gularity at the origin for I — 2, although somewhat complicated in appearance, are equivalent to our 



conditions ( 107 ) if one takes into account Seidel's' definition of A and /i. 



APPENDIX D: STATIC BACKGROUND, AND COMPARISON WITH THORNE ET AL. 

As previous work by Thorne and coworkers, and other authors, was restricted to finding pulsation modes of a static 
background solution, we should point out here how their results fit into our wider framework. 

On a static background the fluid and radial bases coincide, and therefore the comoving and radial coordinates 
coincide. We choose a coordinate system of the form —N'^dt^ + A^dr^ + r^dfl^. (Thorne uses N'^ ~ e'^ and A^ = e^.) 
The background scalars are: 

1 / TTl \ 

U = p = 0, W^—-, iy = Ar{ — +ATTp]. (Dl) 

rA Vr'^ / 

i' and all other time derivatives also vanish. 

Our axial gauge-invariant perturbations, when restricted to RW gauge, are related to the RW gauge axial pertur- 
bations of Thorne by 

u^kA = ^ho, n'^kA^jhi, /? = Ic/t,*, (D2) 

where we have introduced the subscript T to distinguish Thome's variable from our U background scalar. 

We first consider the / > 2 axial perturbations. On a static background we have /3 = 0, and this agrees with 
equation Ut^u = in Appendix B of [|l|. Ref. Q| also states that axial-parity gravitational waves do not couple to 
the stellar matter, because the only matter perturbations that have wave-like solutions (namely sound waves) are 
polar. This is true only on a static background: Eq. ( [7l|) shows that the axial gravitational waves can couple to a 
non- vanishing /3 times an oscillating background coefhcient. 

For the axial I = 1 perturbations, we still have P — 0, and (3 is the only physical perturbation, corresponding to 
the "17" of 1). 

We now consider the I > 2 polar perturbations. For a static background, fi — U — 0, and therefore tp and tp' are 



not present in the wave equations for x and k [see Eqs. (|Al[)-(A2)]. Furthermore, as S^ does not contain ip or ip' , ip 
can be obtained simply by integrating —tp = S^ = 2v{x + k)+x' ■ In this sense, there are only 4 dynamical degrees of 
freedom (not counting the entropy perturbation), while ip plays a passive role. Nevertheless, ■0 can be specified freely 
in the initial data, and in this sense there are really 5 degrees of freedom (not counting the entropy perturbation), 
even on the static background. 

Let us see how this fits into Thome's formalism. Our metric gauge-invariant metric perturbations, when restricted 
to RW gauge, are related to the RW gauge metric perturbations of Thorne by 
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Our matter perturbations a,7, w are related to those of Thorne, Vt, Wt, by 

1 



ip I 1 / P 



^2^ . ^T,. + W't 



p + pj 




r2 2 J 



(D4) 



Again, we have added the subscript T to distinguish two of Thome's variables which might be confused with our 
variables. We have re-obtained all equations (8a-d) ,(9a-c), (C3-4), (C6-7) in ||l[ [except for a wrong r^ factor in (9a)]. 
It is implicit in these equations that Vt, Vt,*, Wt, Wt,*, K and Kt can all be set freely on an initial Cauchy surface. 
Note however that Vt and Wt (undifferentiated) only enter into the equations through their combination in w, and 
we consider only this combination as physical. Thus we count 5 polar perturbation degrees of freedom. 

In discussing the degrees of freedom in more detail we first neglect entropy. The fundamental fluid perturbation 
variable in the framework of Thorne and coworkers is a 3-dimensional displacement vector parameterized by Ut 
(axial) and Vt and Wt (polar). Its time derivative constitutes the fluid velocity perturbation, corresponding to our 
variables /?, a and 7. The displacement generates a perturbation in the particle-number density (and hence in the 
energy density and pressure). 

In our framework, we take the point of view that apart from the velocity perturbations, only the density perturbation 
u! (a linear combination of Vt and Wt) is observable, but not the fluid displacements Ut, Vt and Wt separately. In 
other words, two perfect fluid initial data sets in which the fluid particles are arranged differently, but with the same 
(Eulerian) density and velocity, are considered experimentally indistinguishable. 

The framework of Thorne and coworkers has therefore two degrees of freedom more than ours. They are Ut (axial), 
and a linear combination of Vt and Wt (polar) that is linearly independent of w (for example just Vt). These degrees 
of freedom play only a passive role dynamically, and would be considered "gauge" in our framework. In Thome's 
approach they are considered physical, with the assumption that there exist other properties that provide physical 
"labels" on the particles. 

Finally, we come back to entropy. In our formalism we explicitly add an Eulerian perturbation a to the degrees 
of freedom to account for any kind of entropy perturbation. If the entropy is not constant on the background 
solution, a displacement of the particles would generically produce a (Laplacian) entropy perturbation. The two 
linear combination of Vt and Wt could then be interpreted as the density and entropy perturbations, and both would 
be physical. (But entropy is not considered in [||-g|, and so this is our speculation on an extension of that framework.) 

For the polar I = 1 perturbations, we count 3 degrees of freedom, namely a, uj and 7, while H implies that 4 
functions can be specified freely in the initial data, namely Vt, Vt,*, Wt and Wt.i- This discrepancy is the same as 
in the general I > 2 case. 

Ipser and Price |q| have made a Fourier ansatz for the I > 2 polar perturbations of a static star, with a barotropic 
equation of state p = p{p) for both the background and perturbations, thus addressing the same problem as Thorne 
and coworkers. They find a system of two coupled second-order ODEs for the variables Hq and K of (p3|). This is 
similar to the equations we have obtained for x and fc, minus the equation for ip (which plays only a passive role on 
the static background.) 

APPENDIX E: CHANGE OF FRAME 

Given the fluid basis {uA,nA}, any other orthonormal basis {uA,nA} can be described as a hyperbolic rotation 
(boost) through an angle S,{x^) with respect to the fluid basis: 

ua\(cosH sinhA^MX /«^\ 

UA J \^smh4 cosh4y \^A J \^A J 

This rotation R{£,) induces a transformation in those scalar-variables which are basis-dependent. First, if we 
decompose kAB in the same way as Eq. ( [75[ ) but now using the new basis, we have the same 77, but (j> and ip change 
as 



^;-^(2e)^;) (E2) 



Note that the perturbations a,P,j,LU are basis-independent by definition. 
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The new frame derivatives acting on a scalar / are related to the old derivatives by a single rotation: 

We can change to a new set of projections of v on the new basis and new Ricci-rotation coefficients ji, v: 

Of particular interest besides the fluid basis is the radial basis, where n"^ = v^ /v and u"^ — —e^^fiA- This basis is 
described by ^ = arcth(y), that is, both bases are related by the boost of velocity V. The equations of motion for 
spherical perturbations take their simplest form in the radial base, not the fluid base. In particular, it is useful to 
impose the gauge ip = Q. 

APPENDIX F: ASPECTS OF A NUMERICAL IMPLEMENTATION 



In a previous application of the framework given here ]11| , |12| we have used a numerical evolution scheme that 
explicitly uses the characteristic speeds and is second order in space and time. The set of free evolution equations is 
of the form 

^t+A{r,tf^^=B{r^t)u (Fl) 

in a time coordinate t and a radial coordinate r, where A and B are matrices of rank 3 for the axial perturbations 
and of rank 6 for the polar perturbations. Let V be the matrix of (column) eigenvectors of A. Let A be the diagonal 
matrix composed of the corresponding eigenvalues. Then A = VKV'"^. Let A+ be A with zeros in the place of 
the negative eigenvalues, and A_ be A with zeros in the place of the positive eigenvalues. Let A+ — VK^V^^, 
and let A- — Vh.-V~^, so that A — A+ + yl_. A is quite sparse, and its eigenvalues and eigenvectors can be 
calculated in closed form. With this splitting of A given, we use a finite differencing scheme that is a type of iterated 
Crank-Nicholson scheme, while using the information about left-moving and right-moving modes. We iterate the 
equation 



w^"+i ^ u"; + At 



^"^-'i 2A;^ ^^+ ^^- 2A'r + ■'■ 



(F2) 



Here the auxiliary variable u"+2 is estimated as m" for the first, forward-in-time, step, and then corrected on subse- 
quent steps as u"+2 = (w" -I- m"+-'^)/2. 

The boundary r = does not require special treatment, as u{—r) = ±u(r) for all u, so that ghost grid points with 
r < are available for taking derivatives. The one-sided differencing operators we use here do not give exactly zero 
at r = even if analytically du/dr{0) — 0, but that is consistent: all terms in the finite difference equations combine 
so that numerically, too, m(0) = at all times if u(0) = initially. This also ensures that source terms of the form 
u/r in the evolution equations are well behaved numerically, and do not pose any stability problems. 
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